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ABSTRACT 

We present a parameter that measures the structure of gamma-ray burst time profiles. This parameter is 
based on the statistics of runs and is a good measure of time profile variability. It is shown to be independent 
of burst duration and less sensitive to burst distance than algorithms that depend directly on the intensity of a 
burst. 

Subject heading: gamma rays: bursts 

1. INTRODUCTION 

Gamma-ray burst (GRB) time profiles continue to resist efforts of classification and analysis. As shown by Hurley (1991), these 
profiles show variability on timescales from milliseconds to hundreds of seconds. It is not uncommon in a GRB profile to find peaks 
riding on other peaks or quiescent periods between peaks where the count rates return completely to the level of the preburst 
background. 

If we can parameterize the observed structure, we may be able to classify bursts or find correlations with other parameters. This in 
turn could help place constraints on models of the underlying physical phenomena. Of course, there is no one definition of structure. 
Recently Lamb et al. (1993) used the ratio of the maximum count rates in 64 ms bins to the maximum in the 1024 ms bins as a type of 
burst variability. This definition has the attractive characteristic of applying equally as well to very short bursts as to very long ones. 
On the other hand it is lacking in that, for any one burst, it does not measure the variability over the complete time history. 

In a previous paper (Lestrade et al. 1991), we first presented the idea that a reasonable parameter to measure the variability of a 
profile is the number of times the profile shows long runs of monotonically increasing (or decreasing) bin-to-bin differences. To the 
eye, it is the presence of these monotonic runs that is the signature of what we qualitatively call a highly structured burst. The 
monotonic runs in smooth profiles are short because they are broken up by noise. The statistics of these “runs up,” as they are 
called, although new to the area of GRBs, were heavily studied earlier this century by researchers in agronomy (Wallis 8c Moore 
1941) and meteorology (Besson 1920). This algorithm has the attractive quality of measuring the structure over the entire history. It 
has a shortcoming in that for a given data time resolution, there is a duration limit below which there are not enough bins to count 
runs statistically. 

In a recent paper (Lestrade et al. 1994), we applied this statistic to 173 long-duration BATSE GRB profiles. We found no 
unexpected correlations between structure and other, more easily measurable quantities such as burst angular position, hardness, or 
intensity. Nor did we find any groupings that would indicate more than one source population. In this Letter , we document the 
algorithm used and show its lack of dependence on burst duration and burst V/V t nax . 
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2. RUNS UP AND DOWN 


Consider a gamma-ray burst time profile containing n bins of data. Let the number of counts in the ith bin be 6;. Consider the 
sequence of signs whose ith element is the sign of b i+l — b h (i = 1, 2, ..., n — 1). A subsequence of p consecutive signs (either + or 
— ) is called a run up ( + ) or down ( — ). The notation that we use differs slightly from that used by Knuth (1981) in that, in this Letter , 
a run of length p involves p 4* 1 original bins of data. For example, consider the series of gamma-ray counts registered in the 
following eight 64 ms bins: 501, 502, 504, 510, 503, 502, 503, 501. This series contains four runs (up, down, up, down) of lengths 3, 2, 
1, and 1, respectively. 

Let r p be the number of runs of length p, r’ p be the number of runs of length p or greater, and r be the total number of runs in a 
series. Obviously r\ equals r. Bienayme (1874) was the first to show that the average number of runs in a random series of length n is 
<r> = (2 n — l)/3. The expected values for the number of runs of other lengths are given by Wallis 8c Moore (1941) as 


<Jp> = 2 n 


p 2 + 3p + 1 _ p 3 + 3p 2 - p - 4 
(p + 3)! (p + 3)! 
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Wallis & Moore show that a series of data may be tested for randomness by counting the number and lengths of runs. They 
construct a statistic through the usual procedure of squaring the difference between observed and expected, dividing by the variance 
and summing these ratios, viz.. 
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THe last term of the sum is an integral term over all runs greater than length q — 1. However, as they point out, since these runs are 
not independent, the resultant sums are not distributed according to the usual x 2 statistics. In order to calculate this statistic, we 
need the variances a\r p ) and <T 2 (r' q ). The formulae for these and other covariances have been provided by Levene & Wofowitz (1944). 
They all have the form of <r 2 = an + b where n is the length of the series. Since a and b are roughly the same magnitude and n is large 
for our burst profiles, the an term is more important than b. Equations (4) and (5) give the forms for <j 2 (r p ) and cr 2 (r p ), respectively, 
while Table 1 presents the numerical values of a and b for all values of p up to 7 [the seventh values are for the integral term, <x 2 (r' 7 )] : 
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It might be thought that the number of runs in a series is dependent on the underlying distribution, for example, uniform or 
Poissonian. Wallis & Moore show this not to be the case. While their arguments were based on nonrepeating counts, repetition 
occurs infrequently enough in our profiles so as to not change their conclusion. Certainly the sizes of the bin-to-bin differences are 
affected by the statistical distribution, but the signs of the differences (and consequently the number of runs) are not. Our own 
calculations with Monte Carlo simulations bear this out. 


3. THE PROFILE STRUCTURE PARAMETER, 

Figure 1 presents the number of runs observed in a 60 s interval for (1) two BATSE burst profiles (diamond from burst 2B 910425 
at 2265 s and open circle from burst 2B 920627 at 46,956 s; the latter profile can be seen in Lestrade et al. 1992), (2) BATSE postburst 
background data (open square), and (3) the expected numbers from equations (1) and (2) (line/dot/error bars). The error bars are the 
standard deviations from Table 1. The error bars on the two burst curves are smaller than the diamond and open circle symbols. As 
shown, the numbers seen in the background agree with the statistical predictions. In a burst profile, however, we see increased 
numbers of runs for longer lengths (p t> 3) and usually a smaller number of short runs. It is this difference that we exploit to measure 
a profile’s structure. Following the lead given by equation (3), we have used the expected values and variances to construct a statistic 
that measures the variability of GRB profiles. The statistic is given by 


c _ V ~ <G » 2 , M ~ < r ^» 2 
P + <y\r p ) <rV 7 ) 


(6) 


Note that the r' 7 term accounts for all runs of length 7 or greater. Since we are not strictly trying to measure the deviation from 
statistical chance, but rather trying to measure the structure of a profile by its frequency of long runs, we sum only over terms larger 
than p = 3, and, in addition, we exclude any terms in equation (5) where the observed value is less than the expected value. Not 
doing so would overinflate the variability parameter when a profile showed less structure than random noise at some value of p. 
A highly structured burst will show more long runs (e.g., r 6 , r 7 , etc.) than a smooth one leading to a higher value of S p . Furthermore, 
as will be seen below, S p is insensitive to the duration of the burst. Graphically the parameter S p measures how high the burst curve 
is in Figure 1 above the expected curve for all run lengths greater than p = 3. 

In not counting the shortest runs, our algorithm does miss some short tall “ spikes ” (i.e., runs that are short in duration but tall in 
intensity). However, this is not important for two reasons. Figure 2 is a plot of the number of long runs (p > 3) versus the number of 
short runs but only for runs whose heights are greater than 5.5 a. The first thing to notice is that of the 173 profiles only about 40 
have more than a single tall spike. The other 133 do not have short (tall) runs and therefore do not appear on the graph. (Because we 


TABLE 1 

Coefficients for Variances a 2 ~ an + b 


P 


a 


b 


1 4.236 x tCT 1 -4.819 x 1CT 1 

2 1.127 x l(T l -1.628 x 10" 1 

3 4.192 x 10~ 2 -9.691 x 10 -2 

4 1.077 x 10~ 2 -3.778 x 10" 2 

5 2.004 x 10~ 3 -9.289 x 10~ 3 

6 3.023 x 10' 4 -1.724 x 10‘ 3 

T 4.407 x 10" 5 -3.029 x 10 -4 
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p JG j — Number of runs oflength p in (1) two burst profiles { open circles and diamonds), { 2) postburst background ( open squares), and (3) the expected number in a 
Poissonian background {line/closed circles! error bars) (see text). 

Fig 2. — Number of long runs ( p > 3) vs. number of short runs ( p <, 3) for runs whose heights exceed 5.5 a 


count runs only within the range of 7 90 , we omit the triggering run which is often greater than 5.5 a.) Second, this figure shows that 
of the 40 profiles, there is a good correlation between the number of short runs and the number of long ones. Therefore, statistically 
speaking, the few profiles whose short runs are not counted are also the ones that have the most number of long runs, and so the 
validity of the parameter is preserved. 

4. THE PEAK HEIGHT PARAMETER, H ? 

As mentioned above, one of the appealing attributes of S p is that it is not directly dependent on the heights of spikes, that is, the 
distance to the burster and the intrinsic luminosity function. As a comparison we have constructed a second parameter that is 
directly dependent on peak heights. Following the lead given by equation (5), we formulate a height parameter given by 

H, = £ (h p - </i p » + ( H' p - <//',» , (7) 

P= 1 

where h is the sum of heights for all runs of length p in a particular burst profile, <Ji„> is the height sum observed in a random 
Poissonian profile, and as above, the sum includes only those values where the observed burst height exceeds the noise profile 
height. In short, the value of H p is the sum of all spike (or run) heights above that which you would find in a Poissonian profile of the 

same length. , 

A single-spike burst (Fast Rise, Exponential-Decay or FRED) with no noise would be composed of two runs (one up and one 
down). The value of H would then be twice the height of the FRED. Another burst with two equal FREDs would have twice the 
value of a single FRED ? Because of this definition of H p , longer bursts will have higher values as will be seen below. 

5. CORRELATION WITH BURST DISTANCE AND DURATION 

Figure 3 shows that both S. and H depend, to different degrees, on the burst distance or V/V mtl . The structure parameter, being 
dependent on the number of runs and not their heights, shows slightly less correlation with distance than does the height parameter. 
The profiles appearing in these correlations are the same 173 long-duration bursts from Lestrade et al. (1992). 
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Fig. 3.— Burst V fV m „ vs. S, and H 0 
Fig. 4. — Burst duration vs. S p and H p for normalized profiles (see text) 
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Fn order to remover the dependence on distance, we normalized the burst profiles by reducing all of them to the same value of 
C max . The algorithm to handle the reduction and partial addition of Poissonian noise was obtained from Norris & Nemiroff (1994). 
Figure 4 presents the graph of duration versus S p and H p for these normalized profiles. As a measure of duration, we use the T 90 
values from Fishman et al. (1994). With the intensities and noise biases rendered uniform, we see very clearly the dependence of H p 
on duration. The structure S p , however, does not show a correlation with burst duration. 

6. CONCLUSION 

The parameter S p is a measure of the variability of a burst profile. It has several properties that make it attractive as a statistic. In 
addition to being easily formulated and straightforward to. calculate, it is independent of background count rate, trigger integration 
time, and burst duration, and can be made independent of peak intensity. Furthermore, it is less dependent on burst distance than 
are parameters that depend directly on the heights of profile features. 

The author wishes to thank Brad Shaefer, David Palmer, and Jay Norris at GSFC for several helpful comments. This research 
was supported by NASA grant SUB 92-194. 
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